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Abstract It's well-known that there is a so-called high-level error bound 
for multiquadric and inverse multiquadric interpolations, which was put for- 
ward by Madych and Nelson in 1992. It's of the form \f(x) — s(x)\ < A^ \\f\\ h 
where < A < 1 is a constant, d is the fill distance which roughly speaking 
measures the spacing of the data points, s(x) is the interpolating function 
of f(x), and h denotes the multiquadric or inverse multiquadric. The error 
bound converges very fast as d — > 0. The constant A is very sensitive. A slight 
change of it will result in a huge change of the error bound. Unfortunately 
A can not be calculated, or even approximated. This is a famous question in 
the theory of radial basis functions. The purpose of this paper is to answer 
the question. 
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1 Introduction 

Let h be a continuous function on M. n which is conditionally positive def- 
inite of order m. Given data (xj, fj) , j — 1, N, where X = {xi, x^ } 
is a subset of points in W 1 and the fj s are real or complex numbers, the 
so-called h spline interpolant of these data is the function s defined by 



N 



'1.1) s(x) = p(x) + Cjh( 
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where p(x) is a polynomial in V m -\ and the c/s are chosen so that 



N 




for all polynomials q in V m -\ and 



N 



(1.3) 



Xi Xj) fi, i !)•••) N. 



Here V m -\ denotes the class of those polynomials of 1R™ of degree < m — 1. 

It is well known that the system of equations (1.2) and (1.3) has a unique 
solution when X is a determining set for V m -\ and h is strictly conditionally 
positive definite. For more details please see [5]. Thus, in this case, the 
interpolant s(x) is well defined. 

We remind the reader that X is said to be a determining set for V m -\ if 
p is in V m -i and p vanishes on X imply that p is identically zero. 

In this paper h is defined by formula 



where \x\ is the Euclidean norm of x, T is the classical gamma function 
and P, c are constants. The function h is called multiquadric or inverse 
multiquadric, respectively, depending on (3 > 0, or (3 < 0. 

In [6] Madych and Nelson obtained bounds on the pointwise difference 
between a function / and the h spline which agrees with / on a subset X 
of R n . These estimates involve a parameter that measures the spacing of 
the points in X and are O (d £ ) as d — > where I depends on h. Later in 
[7] they found that for multiquadrics and inverse multiquadrics, the estimate 

can be improved to O as <i ^ 0, where A is a constant which satisfies 

< A < 1. The conditions on / are the same as those in [6]. 

1.1 A Bound for Multivariate Polynomials 

A key ingredient in the development of our estimates is the following 
lemma which gives a bound on the size of a polynomial on a cube in IR ra in 
terms of its values on a discrete subset which is scattered in a sufficiently 
uniform manner. We cite it directly from [7] and omit its proof. 

Lemma 1.1 Forn = 1,2, define j n by the formulas 71 = 2 and, if n > 1, 
7„ = 2n(l +7«-i). Let Q be a cube in M, n that is subdivided into q n identical 



(1.4) 



h(x) :=r(-^)(c 2 + |x| 2 )f, (3eR\2N: 
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subcubes. Let Y be a set of q n points obtained by selecting a point from each 
of those subcubes. If q > j n (k + 1), then for all p in Vk 

sup \p{x)\ < e 2n 7n(fc+i) sup \ p ( y )\ _ 

xeQ yeY 

1.2 A Variational Framework for Interpolation 

The precise statement of our estimates concerning h splines requires a 
certain amount of technical notation and terminology which is identical to 
that used in [6] . For the convenience of the reader we recall several basic 
notions. 

The space of complex-valued functions on M, n that are compactly sup- 
ported and infinitely differentiable is denoted by V. The Fourier transform 
of a function in T> is 

0(0 = J e- t<x ^><f)(x)dx. 
A continuous function h is conditionally positive definite of order m if 

h(x)(j)(x) * (f>(x)dx > 

holds whenever = p(D)if) with ip in D and p(D) a linear homogeneous 
constant coefficient differential operator of order m. Here 4>{x) = (f)(—x) and 
★ denotes the convolution product 

4>1 * 02 (t) = J 4>l{x)4>2{t - X)dx. 

As pointed out in [6] , this definition of conditional positive definiteness is 
equivalent to that of [5] which is generally used. 

If h is a continuous conditionally positive definite function of order m, 
the Fourier transform of h uniquely determines a positive Borel measure fi 
on R n \ {0} and constants a r ,\r\ =2m as follows: For all vp G V 

(i.5) J h(x)^x)dx = J | m - m E d ^(°)t7 I mo 



|7|<2m 



|7|<2m 
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where for every choice of complex numbers c a ,\a\ — m, 

a a+/3 c«c^ > o. 

|o|=m |/3|=m 

Here x is a function in V such that 1 — has a zero of order 2m + 1 at 
£ = 0; both of the integral J 0<|{|<1 |£| 2m are finite. The 

choice of x affects the value of the coefficients a 7 for | --y J < 2m. 

Our variational framework for interpolation is supplied by a space we 
denote by C h , m - If 

V m = |(/> G V : y x a (j)(x)dx = /or a// |a| < m| , 

then C^ m is the class of those continuous functions / which satisfy 

(1.6) 



f(x)4>(x)dx 



< c (f) W K x ~ y) ( P(x)<P(y)dxdy 



for some constant c(/) and all in P m . If / e C h ^ m , let denote the 

smallest constant c(f) for which (1.6) is true. Recall that ||/||^ is a semi- 
norm and Ch, m is a semi-Hilbert space; in the case m = it is a norm and a 
Hilbert space respectively. 

2 Main Results 

We first recall that the function h defined in (1.4) is conditionally positive 
definite of order m = if (3 < 0, and m — [|] if (3 > 0. This can be found in 
[8] and many relevant papers. Then we have the following lemma. 

Lemma 2.1 Let h be as in (1.4) and m be its order of conditional positive 
definiteness. There exists a positive constant p such that 

(2.1) / 1^1" dp (0 < (V2) n+P+1 ■ (vVf +1 • na n ■ c^ k ■ A • p k ■ k\ 

for all integer k > 2m + 2 where p is defined in (1.5) , a n denotes the volume 
of the unit ball in R n , c is as in (1.4) , and A is a positive constant. 
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Proof. Let K v denote a modified Bessel function of the second kind. Then 

n+0 

n + 

= 2 ^{jc) 2 ■ J^t^ -K^{c\i\)di 

n + 

~ -V-27T2- — / fr 2 ■ - di 

^2 \2cJ J R J q] v^Wl-e^ 



n + 

"°° j_a r "-l 

j < ._n±0 T 



K ■ 27T2 • ( j ■ n - a n I r K ■ dr 



V2 ~ V 2c / Jo yfc\r\-e c \ r \ 



■^7= • 27T2 • (2c) 2 • — • n • a n dr 

v^F ft n , n ,n±A 1 1 f°° r k+^ 



n-0-3 
2 



= -^-27T2-(2c) 2 • — = • n • a n • ^— r f dr 

V2 K ' ^fc c k +^Jo e r 

k 

= 2 2 . tt^ - • n • 0K n • c^ - ^ / — dr where k = k + — — — — -. 

J e r 2 

Note that if (3 < 0, then m = and k > 2m + 2 = 2. This implies k' > 0. If 
P > 0, then m = ff ] and fc > 2m + 2 = 2 [§] + 2. This implies k' > 0. In 
any case k' > 0. 

Now we divide the proof into three cases. Let fc" = \k'~\ which is the 
smallest integer greater than or equal to k' . 

Casel. Assume k" > k. Let k" = k + s. Then 

/ // 

/oo fc /-oo k 

—dr < J —dr = k"\ = (k + s)(k + s — 1) • • • (k + l)k\ 

and 

poo k'+l roc k"+l 

/ -^dr< —^dr = (k" + 1)1 = (k + s + l)(k + s)--(k + 2)(k + l)kl 
Jo eV Jo eT 

Note that 

(k + s + 1)0 + s) ■ ■ ■ (k + 2) _ k + s + 1 
(k + s)(k + s-l)---(k + l) ~ k+1 
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(i) Assume (3 < 0. Then m = and k > 2. Tm's gives 

A; + s + 1 < 3 + s 



fc+ 1 



3 

oo „fe"+l 



/•oo fe +1 

/ -^rfr< A -p fe+1 -(A; + l)!. 



i/ J °° ^dr < A • p k • k\. The smallest k" is k^ — 2 + s. iVo-u;, 

/ —dr = k \ = (2 + s)(2 + s- 1) • • • (3) • k\ where k = 2 

= ( 2 + '" 2 + 7 1 >-' 3 >.AK* = 2) 



= A • p 2 ■ 2! where A = 



(2 + s)(2 + s- 1) • • • (3) 



P 2 



.ft follows that J °° ^rdr < A • p fc • k\. for all k > 2. 

(ii)^lssume /? > 0, Then m= [~f ] and > 2m + 2. T/i?s gives 



k+s+l < 2m+3+s_^ s 



k + 1 ~ 2m + 3 2m + 3 



I 



oo r .k +1 

dr< A • p fc+1 • (k + 1)! 



i/ J °° < A ■ p k ■ k\. The smallest k" is k" Q =2m+2+s when k=2m+2. 

Now, 

/o°° 4dr 

= fcj! = (2m + 2 + s)(2m + 1 + s) • • • (2m + 3)(2m + 2)! 
( 2m + 2 + S )(2m+l + g )---(2m + 3) 2m+2 

_ p 2m+2 "P "^m + Zj. 

2m+2 nM , a (2m + 2 + s) (2m + 1 + s) • • • (2m + 3) 
= A • p 2m+2 • (2m + 2)! where A = - ^ 2m+2 
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It follows that J °° r —dr < A p k k\ for all k>2m + 2. 

Case2. Assume k" < k. Let k" = k — s where s > 0. Then 
/ // 

roo r k roo ^fc ^ 

/ — dr < / — dr = k"\ = (k - s)\ = — : - • k\ 

Jo e r -J e y 1 fc(fc-l)...(fc-s + l) 



and 

r .k +1 poo r k + 1 

dr 



/ —dr < 

Jo e r J 



& 



= '*' + 1 " = ( t - +1 >' = ( t+ l) t ...( t -. + 2) -( t+1 >'- 

Note that 

1 / 1 | fc(]fe - 1) • • • (k - s + 1) k-s+1 



(k + l)k ■ ■ ■ (k - s + 2) / k(k - 1) • • • (k - s + 1) J (jfc + l)fc • • • (fc - s + 2) fc + 1 

(i) Assume (5 < 0, £/ien m = ane? k > 2. Since k" = k — s > 1 /ioWs /or 
all k > 2, it must be that s — 1. TTras 

k o _l_ 1 O 1 

— — = 1 - = 1 - < 1 for all k> 2. 

k+1 k+1 k + 1 ~ 

Let p — 1. TTien 

"oo „fe +1 /-oo k 



/ ^rfr< A p fe+1 (A: + 1)! i/ / "—dr < A p k k\. 
Jo e Jo e 

The smallest k" is k' ' = k — s = 2 — s. Now, 
f°° r k o 

/ —dr = k \ = (2-s)\ = l\ = 1 
Jo e 



= -k\ where k = 2 
2 

= A p k k\ where A = -. 

It follows that f Q °° r —dr < A p k k\ for all k>2. 

iii) Assume f3 > Then m=\^ ] and k > 2m + 2. This gives 

k ~ S + 1 = 1 g— < 1. 

fc+ 1 ifc + 1 _ 
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Let p = 1 Then ^-^dr < A • p k+1 ■ (k + 1)! if f °° r -^dr <A -p k -k\. 



The smallest k is ko = 2m+2. Hence the smallest k" is k' ' =ko—s = 2m+2—s. 
Now, 

f°° r k o 

/ —dr = k \ = (2m + 2-s) = (k - s)\ 
Jo eT 

' (ko)\ 



k (k - 1) • • • (k - s + 1) 
A • p ko k \ where A = 



(2m + 2)(2m + 1) • • • (2m - s + 3) 
It follows that / °° r —dr < A p k k\ for all k>2m + 2. 
Case3. Assume k" = k. Then 

roc r k' roo r k" roo ^ +1 

/ — dr < / — dr = k\ and / ^-dr < (k + 1)!. 
Jo eT Jo eT Jo eT 

Let p = 1. Then 

roo r k 

/ — dr < A p k k\ for all k where A = 1. 
Jo e r 

The lemma is now an immediate result of the three cases. □ 

Remark. For the convenience of the reader we should express the constants 
A and p in a clear form. It's easily shown that 

(a) k" > k if and only if n — (3 > 3; 

(b) k" < k if and only if n — (3 < 1; 

(c) k" = k if and only if 1 < n — /? < 3 

where fc and A; are as in the proof of the lemma. 

We thus have the following situations, 
(a) n-(3 > 3. Let s = [^=f^] . Then 

(ii)i/ /?>0,p=l + ^ i and A = (^2+^(2^+1+*)... ( 2m+3 ) 
where m = [|] . 
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(b) n-/3<l. Let s = - \^=^] ■ Then 



1. 

2' 



(i)if (3 < 0,p = 1 and A 

P > 0,p = 1 and A = ( 2Tra+2)(2m+ l ) ,„ (2m _ s+3 ) ; w/»ere 



(c) 1 < n-/3 < 3. Then p = 1 and A = 1. 

Before introducing our main theorem, we need the following lemma which 
is cited directly form [7]. 

Lemma 2.2 Lei Q, 1", and j n be as in Lemmal.l. Then, given a point x in 
Q ,there is a measure a supported on Y such that 

I P(y)da(y) = p(x) 



for all p inVk, and 



[ d\a\ (y) < e 2n ^ k+1 \ 



Now we need another lemma. 
Lemma 2.3 For any positive integer k, 



k\ - ■ 

Proof. This inequality holds for k — 1 obviously. We proceed by induction. 



y/[2(fc + l)]! _ ^(2^ + 2)! _ y^2fc)! y/(2fc + 2) ■ (2fc + l) 
(k + l)\ ~ k\(k + l) ~ k\ k + l 

<} /m VW+W <2k (2k + 2) 
~ k\ k+l ~ k + l 

□ 

Because of the local nature of the result, we first restrict our attention to 
the case where x lies in a cube. 

Theorem 2.4. Suppose h is defined as in (1.4) and m is its order of con- 
ditional positive definiteness. Let p be its corresponding measure as in (1.5) 
.Then, given a positive number bo, there are positive constants 5q and A,0 < 
A < 1, which depend on 6 for which the following is true: 
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If f G Ch, m an d s is the h spline that interpolates f on a subset X of R n , 
then 

(2.2) |/(x) - s(x)\ < 2^ • vr^ ■ ^cT n • c f • • A* • ||/|| fc 

holds for all x in a cube E provided that (i)E has side b and b > &o>(ii)0 < 
5 < 5o and (iii) every subcube E of side 5 contains a point of X. Here, a n 
denotes the volume of the unit ball in R n and c, A are as in (2.1) . 
The numbers 5 and X can be expressed specifically as 



2 \ 3C 7 n 



where 



6 ° 3C 7n (m + l)' A= V3 

C = max {2p'^ie 2n ^, ^-X , p = -. 

{ 36 J c 



The number p can be found in the remark immediately following Lemma2.1, 
and 7 n is defined in Lemmal.l. 

Proof. First, let p, 7„, and 6 be the constants appearing in (2.1) , Lemmal.l. 
and Theorem2.4, respectively. Let 



B = 2p \fne 2n '~ 1n and C = max < B, — — 1 where p 

{ 36 J 



P 

c 



Let 

&= 1 



3C ln (m + 1) ' 

where m is the order of conditional positive definiteness of h. 

Now, let x be any point of the cube E and recall that Theorem4.2 of [6] 
implies that 

(2-3) |/(:r) -*(*)!<<* II/IL / \y-x\ k d\a\(y) 

whenever k > m, where a is any measure supported on X such that 

(2.4) / p(y)da{y) = p(x) 

for all polynomials p in Vk-i- Here 

1 

whenever k > m. By (2.1) , for all 2k > 2m + 2, 
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(2-5) c k = {j Rn ^d^)\ 



1 n+p + 1 n + 1 , 0-2k /— ■ n — — 



n + /3 + l n + 1 £ . 

<2 4 . 7T 4 . C 2 • c fe • a/Ao • ( 2 P) 

due to Lemma2.3. 

To obtain the desired bound on \f(x) — s(x) \ , it suffices to find a suitable 
bound for 



I = Ck / \y- x\ k d\cr\ (y). 



This is done by choosing the measure a appropriately. We proceed as follows: 
Let 5 be a parameter as in the statement of the theorem. Since 5 < S 
and < 3Cr n 5 < ^L_, We may choose an integer k > m + 1 so that 

1 < 3C~f n k5 < 2. 

Note that 'jnkd < bo for such a k. Let Q be any cube which contains x, has 
side jnkd, and is contained in E. Subdivide Q into ( n f n k) n congruent subcubes 
of side S. Since each of these subcubes must contain a point of X, select a 
point of X from each such subcube and call the resulting discrete set Y. By 
virtue of Lemma2.2 we may conclude that there is a measure a supported 
on Y which satisfies (2.4) and enjoys the estimate 



(2.6) / d\a\{y)<e 



2rfy n k 



We use this measure in (2.3) to obtain an estimate on I. 

Using (2.5) , (2.6) , and the fact that support of a is contained in Q whose 
diameter is y/n^ n k5, we may write 



(2.7) / < 2^r~ • vr^ • v /^ • • c~ k • • (2p) k (^i ln k5) k e 



< (C ln kS) k (2^ ■ tt^ • ■ J ■ ^A~ ). 

Since 

2 1 

Cj n kS < - and k > 



k p 2n^ n k 



3 " 3C 7n <T 
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(2.7) implies that 



I < 



i 

- ) 3C 7 . 

3 J 



4 n+P+l n+l §_ 

■ (2 4 . 7T 4 . y/na n ■ c 2 



Hence we may conclude that 

|/(x)- S (x)|<2 I1± l ±i -vr^ 



where 

A 



2 \ 3C 7n 

3. 



□ 

What's noteworthy is that in Theorem2.4 the parameter 5 is not the 
generally used fill distance. For easy use we should transform the theorem 
into a statement described by the fill distance. 

Let 

d (Q, X) = sup inf \y — x\ 

yen x ^ x 

be the fill distance. Observe that every cube of side S contains a ball of radius 
|. Thus the subcube condition in Theorem2.4 is satisfied when 5 = 2d(E, X). 
More generally, we can easily conclude the following: 

Corollary2.5 Suppose h is defined as in (1.4) and m is its order of condi- 
tional positive definiteness. Let /j, be its corresponding measure as in (1.5) . 
Then, given a positive number b , there are positive constants d and A , 
< A < 1, which depend on b for which the following is true: If f G (\ m 
and s is the h spline that interpolates f on a subset X of R n , then 

(2.8) \f{x) - s{x)\ < 2^ ■ ■ ^cT n • J • y/%> • (A')* • ||/|| fc 

holds for all x in a cube E C Jl, where Q is a set which can be expressed 
as the union of rotations and translations of a fixed cube of side bo, provided 
that (i)E has side b > bo, (ii)0 < d < d and (in) every subcube of E of side 
2d contains a point of X. Here, a n denotes the volume of the unit ball in R n 
and c, A are as in (2.1) . Moreover d = y and \' = V\ where 5q and X 
are as in Theorem2.\. 

Proof. Let d = f and 5 = 2d. Then < d < d iff < 5 < 5 . Our 
corollary follows immediately by noting that = = v / A d = (X')^. □ 
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Remark. The space Ch, m probably is unfamiliar to most people. It's in- 
troduced by Madych and Nelson in [5] and [6] . Later Lull made characteriza- 
tions for it in [1] and [2] . Many people think that it's defined by Gelfand and 
Shilov's definition of generalized Fourier transform, and is therefore difficult 
to deal with. This is wrong. In fact, it can be characterized by Schwartz's 
definition of generalized Fourier transform. The situation is not so bad. 
Moreover, many people mistake C h , m to be the closure of Wu and Schaback's 
function space which is defined in [9] . This is also wrong. The two spaces 
have very little connection. Luh also made a clarification for this problem. 
For further details, please see [3] and [4] . 
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